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Abstract 

As long as we neglect backreaction, the Hawking temperature of a given black hole would 
not depend upon the parameters of the particle species we are considering. In the semiclassical 
complex path analysis approach of Hawking radiation, this has been verified by taking scalar and 
Dirac spinors separately for different stationary spacetime metrics. Here we show, in a coordinate 
independent way that, for an arbitrary spacetime with any number of dimensions, the equations 
of motion for a Dirac spinor, a vector, spin-2 and spin-| fields reduce to Klein-Gordon equations 
in the WKB semiclassical limit. We then obtain, under some suitable assumptions, the complex 
solutions of those resulting scalar equations across the Killing horizon of a stationary spacetime 
to get a coordinate independent expression for the emission probability identical for all particle 
species. Finally we consider some explicit examples to demonstrate the validity of that expression. 
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1 Introduction 



The semiclassical tunneling method pQ- [7] is an alternative approach to model particle creation by 
black holes [8] . The basic scheme of this method is to compute the imaginary part of the 'particle' 
action which gives the emission probability from the event horizon. From the expression of the 
emission probability one identifies the temperature of the radiation. The earliest works in this 
context can be found in [TJ [2] . Following these works an approach called the null geodesic method 
was developed O 0] . There exists also another way to model black hole evaporation via tunneling 
called complex path analysis El [7] which we discuss here. This method involves writing down, 
in the semiclassical limit Ti — > a Hamilton- Jacobi equation from the matter equations of motion, 
treating the horizon as a singularity in the complex plane (which is a simple pole for all known 
solutions) and then complex integrating the equation across that singularity to obtain an imaginary 
contribution for the particle action. 

Both of this two alternative approaches have received great attention during last few years. It is 
noteworthy that since both of these methods deal only with the near horizon geometry, they can be 
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very useful alternatives particularly when the spacetime has no well defined asymptotic structure 
or infinities [3]. 

As far as we neglect the backreaction of the matter we are considering, the temperature of the 
radiation or the Hawking temperature should not depend upon the parameters (e.g. mass, spin, 
and charge) of the particle species. The Smarr formula for black hole mechanics predicts that this 
temperature is proportional to the surface gravity of the event horizon for a stationary black hole 
with a Killing horizon. 

The complex path analysis approach has been successfully applied to scalar emissions as well 
as to spinor emissions separately for a wide class of stationary black holes giving the expected 
expressions of Hawking temperatures that were predicted by the Smarr formula. To tackle Dirac 
equation in this approach the usual method has been employed, i.e., finding a proper representation 
of the general 7 matrices in terms of the Minkowskian 7's and the metric functions and then making 
the variable separation. For an exhaustive review and list of references on this see e.g. jlOj . See 
also e.g. |llj-|20) for some recent issues concerning the tunneling approach. 

Thus, the universality of the Hawking temperature has been proved case by case for a wide 
variety of black holes via the complex path method. Can we prove this universality from a more 
general point of view? 

In particular, in this paper we shall show that for the Dirac spinors we do not need to work 
with any particular representation of the 7 matrices in the semiclassical framework. In this work 
we wish to point out, in a coordinate independent way that in any arbitrary spacetime with any 
number of dimensions, the equations of motion for a Dirac spinor, a vector, spin-2 meson and spin-| 
fields reduce to Klein-Gordon equations in the semiclassical limit ft — > for the usual WKB ansatz. 
The equations for a charged Dirac spinor reduce to that of a charged scalar. This clearly shows 
that at the semiclassical level all those different equations of motion of various particle species are 
equivalent and it is sufficient to deal with the scalar equation only. We shall also present, for a 
stationary spacetime with some assumed geometrical properties, a general coordinate independent 
expression for the emission probability and the Hawking temperature which is characterized by the 
black hole parameters itself (Eq. ([23)) ). We further consider some explicit examples to demonstrate 
that our formula indeed gives the expected Hawking temperature in terms of the horizon's surface 
gravity. 

Thus the semiclassical complex path method gives us a way in which we may treat the different 
spin fields in an identical footing, giving the same Hawking temperature and thereby proving the 
universality of the Hawking temperature for stationary black holes from a very general point of 
view. 

The paper is organized as follows. In the next section we shall deal with Dirac spinors (neutral 
and then charged) to show that the equations reduce to that of scalars in the semiclassical limit 
for the WKB ansatz. In Sect. 3, we shall explicitly expand the resultant scalar equation in a 
coordinate independent way in the near horizon limit for a stationary black hole with a Killing 
horizon, and shall present a general expression that gives the emission or absorption probabilities. 
We shall illustrate the validity of this expression by taking a few explicit examples. In Sect. 4, 
we shall also demonstrate that similar results hold also for the vector, massive spin-2 and spin-| 
fields. Finally we shall discuss our results. 

We shall take G = 1 = c, but shall retain ft, throughout. 

2 Reduction of the semiclassical Dirac equation into Klein- 
Gordon equation 

Let us then start by considering a spacetime of dimension n, and a metric g a b defined on it, at least 
in our region of interest. We consider the Dirac equation 
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V a is the spin covariant derivative denned by := (d a + T a ) "J, where T a are the spin connection 
matrices. The matrices "f a (x) are the curved space generalization of the Minkowskian 7W. We 
expand 7" in an orthonormal basis, 7° = 7^e? % : /i = 0, 1, 2, . . . , (n — 1). Also, g 



V 



Here the Greek indices within bracket denote the local Lorentz indices and J7<M(") is the 



inverse metric corresponding to the n-dimensional Minkowski spacetime. The 7" satisfy the well 



known anti-commutation relation: {7 



'} = 2r)^^I, where I denotes the identity matrix. 



The expansion of 7° in terms of the orthonormal basis {e"^}, and the anti-commutation relation 
for 7^'s give 

{ 7 a , 7 6 } = 2g ab l. (2) 
Now we square Eq. ([T]) by acting with ij h Vb on both sides from left, producing 



~ (7V + 7°7 b ) V b V a * + ~ (7 V - 7°7 b ) (V 6 V a - V a V h ) 



*+(7 b V b7 Q )V a * = -^*. (3) 



But the commutator of two covariant derivatives acting on ^ is proportional to the Riemann tensor, 



(7V - 7 a 7") (V h V Q - V Q V b ) 



7 7 



7 b 7" 



bed 



\t. Using this fact and the 



anti-commutation relation for 7° (Eq. @), Eq. ([3]) becomes 



V Q V a * + - [ 7 a , 7 fo ] R abcd [ 7 C , 7 d ] * + ( 7 6 V 67 a ) V a * 



(4) 



We will look at Eq. Q semiclassically. We choose the usual WKB ansatz for a spin-'up' particle 



A(x) 


B(x) 




(•5) 



and substitute into Eq. ([4]). Since we are neglecting backreaction, the components of the Riemann 
tensor are independent of fi. Then it is clear that in the semiclassical limit Ti — > 0, on the left hand 
side only the first term survives because only this one contains some double derivatives of O (h~ 2 ) . 
The single derivative terms coming from the Laplacian will certainly not survive in the semiclassical 
limit (which is true for an actual scalar equation also), but we shall formally keep the Laplacian 
V a V a intact till later when we shall discuss its expansion explicitly. Thus in the semiclassical limit, 
the WKB ansatz ((5j) implies Eq. (|4]) can be effectively represented by two Klein-Gordon equations 
for spin-'up' particles 



V a V a ^ + -5-$ = 0. 

h 



(6) 



Similar result holds for a spin-'down' particle also. 

If we consider a Dirac particle with a charge e coupled to a gauge field A a , the spin covariant 
derivative V„ in Eq. (P) is replaced by the gauge covariant derivative V a = V a — j^A a such that 
the equation of motion becomes 



(7) 



We now apply from the left (i7 b V b 
([3]) we obtain 

V a V a * + i [7", 71 Rabcd [l C - 



ie 
J 



^7 b j4 b ) on both sides of this equation. Using Eq.s ([2]) and 
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p 2 Oj'p 

(7 b V b7 a ) V a * - -^A b A b ^ + --A a V a y 
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(7 b V b7 Q ) A a + - [ 7 a , 7 b 



F„ 



(V a A a ) 



(8) 
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where F ab — \7 a A b — Vf,A a , We now substitute the ansatz (Eq. (JSJ) ) into Eq. © and take the 
semiclassical limit ?i — » 0. We see that in this limit Eq. © can formally be represented by 

V Q V a * - + ^ a V a * + = 0, (9) 

h n h 

each of which effectively has the form of the equation of motion of a charged scalar. 

What have we seen so far? We have dealt with neutral and charged Dirac spinors and have 
explicitly shown in a coordinate independent way that, for the semiclassical WKB ansatz all those 
equations of motion are equivalent to that of scalars in any arbitrary spacetime of dimension n. So 
it is clear that the single particle Hawking radiation will be identical for Dirac spinors and scalars 
for any given black hole. 

We shall also show explicitly in Sect. 4 that similar conclusions hold for Proca, massive spin-2 
and spin- 1 fields. But before that we wish to discuss the explicit expansions and the near horizon 
limits of Eq.s ©, ([9]) in a stationary spacetime containing black hole. We shall address only the 
charged Dirac spinor (or equivalently, charged scalar, Eq. ([9])). The other case will be equivalent 
to setting e = in Eq. ([9]). 



3 Hawking temperature for a stationary black hole with 
Killing horizon 

We wish to present in the following a general coordinate independent expression for the emission 
or absorption probability from a stationary black hole with some assumed geometrical properties. 
Let us first list the definitions and assumptions we make. 

We consider an n-dimensional stationary spacetime containing a black hole with a Killing hori- 
zon T-L. We assume that the spacetime can be foliated into a family of hypersurfaces E, orthogonal 
to a vector field x a ■ The hypersurface is spacelike everywhere except at the horizon (%), which is 
defined to be an (n — 1) dimensional null hypersurface. So, \ a is orthogonal to a null hypersurface 
over Ti and hence x° is itself null over Ti. Everywhere else \ a is timelike. 

Since T-L is a Killing horizon, the vector field x° becomes a null Killing vector field, say Xh j over 
T-L. \ a is not necessarily a Killing field everywhere, but it is Killing at least over Ti 

X a \u = Xh : V (aXHh ) - 0, XnXHa = ~P 2 \n = 0. (10) 
We now write the spacetime metric g a b as 

9ab = -fi~ 2 XaXb + \~ 2 R a Rb + lab, (11) 

where R a is a spacelike vector field orthogonal to x a , an d A 2 is the norm of R a . ^ a b is the non-null 
spacelike portion of the metric perfectly well behaved on or in an infinitesimal neighbourhood of 
the horizon. 

Let us denote the Killing fields of this spacetime by (£„, {fia}), where i — 1,2... to. Let 
£ a be the timelike Killing field and be the spacelike Killing field(s). We assume that the 

hypersurface orthogonal vector field x a (which is orthogonal to {4> l a } and any other spacelike field), 
can be written as a linear combination of all the Killing fields 

Xa=Za + a i {x)4? a , (12) 

where repeated indices are summed over and {a l (x)} are smooth functions. Then, using Killing's 
equation we have V( a Xfc) = ( / , ( a V{,)a l (x). Thus we have 

X a X h ^aXb = -\x a V a P 2 - xV^VbC^z) = 0. (13) 
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Eq. ([T5)) shows that V a /3 2 is everywhere orthogonal to \ a an d hence it is spacelike when x a is 
timelike. So, we may choose R a = V a /3 2 in Eq. (ITT1) . 

To look at the behaviour of V a /3 2 over the horizon, we recall that over the Killing horizon % 
EH [22] 

V a /3 2 = -2 KX Ha, (14) 

where k is a function. Since by definition %h i s nuu hypersurface orthogonal at the horizon, it turns 
out that k is a constant over the horizon [21]. Eq. (|14|) shows that V a /9 2 is null over T~L. However, 
the choice R a — V a /? 2 is not unique, we could have multiplied V a /3 2 by some non-diverging function 
over T-L, even some positive power of /3. But we shall retain this choice for convenience. 
Let R be the parameter along R a . Then using Eq. (fTi|) we have over T-L 

# a V /3 2 = ^ = -4 K 2 /? 2 , (15) 

which implies over T-L 

/3 2 = er^ R . (16) 

With the choice of R a we have made, it is clear that the metric (ITT1) becomes doubly degenerate over 
%. Note that Eq. (fTTj) can readily be realized, in its doubly degenerate form, for a static spherically 
symmetric black hole by employing the usual (t, r*) coordinates, where r* is the Tortoise coordinate. 
We shall be more explicit about R when we shall go into specific examples. 

The assumption of stationarity and Killing horizon would help us to provide a meaningful notion 
of the 'particle' energy [21] . 

For n > 4, the uniqueness and other general properties of black holes are not very well under- 
stood and there may exist more general stationary black holes. However, we shall show below that 
for known stationary exact solutions, those assumptions will be sufficient. 

Let us now expand Eq. © with the decomposition ([TTj). The single derivative terms do not 
contribute in the fi — > limit we are concerned with and the equation explicitly becomes 

A 2 (x a d a I - eff - f3 2 (R a 8 a I + eg) 2 - ((3X) 2 [ lab d a Id h I + e 2 lab A a A h - 2e lab A a d b I + m 2 ] = 0, (17) 

where / = — x a ^4 a , and g — R a A a . Here it is clear that had we multiplied R a by a function h(x) 
non-diverging over %, we would have multiplied Eq. (|17[) only by an over all factor h 2 {x). 

Now we shall look Eq. (|17[) in the near horizon limit. By our assumption the metric functions 
"fab are well behaved over the horizon. So, r y a bA a A b is non divergent over T-L. Also, examples with 
g ^ seem to be unknown in the literature. So, we shall set g = in Eq. (fT7|) and write Eq. (fl7|) 
in the near horizon limit as 

A 2 ( X a d a I - eff - /3 2 (R a d a Tf - ((3X) 2 [ lab d a Id b I - 2e lab A a d b l] = 0. (18) 

To further simplify Eq. ([T5|) . let us choose an orthogonal basis {m 1 *}™' 2 for j ab . Let 9i be the 
parameter along each m". Let us consider the first term within the square brackets. This is basically 
a sum of the squares of (n — 2) Lie derivatives: -^i{£ mi I) 2 + ^(£m 2 I) 2 ■ ■ •> where m 2 is the norm 
of each m". By our definition, those norms are non-zero finite over %. Since / is a scalar those Lie 
derivatives are basically partial derivatives : £ mi I = dg i I. 

We shall now check whether the terms within the square bracket in Eq. (|18|) are divergent over 
H. Let us suppose that close to H, if possible the following divergence occur 

labd *Id b I=^, (19) 
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where D(x) is bounded over or close to T-L and independent of (3 at leading order. Then Eq. (|15[) 
implies that D(x) is also independent of R over T-L 



£ R D(x) = 0. (20) 

Also by our choice R a — V a /3 2 , whose norm is A 2 , vanishes over % as 0((3 2 ) (Eq. So the 

function £>(ie) is also independent of A in the leading order over T-L. Since the metric functions jab 
are well behaved over T-L, the divergence of "f ao d a Id b I arises from the Lie derivatives (dgJ) 2 . For 
simplicity we shall suppose that the divergence comes from a single Lie derivative which is the i-th 
one. We can easily generalize our analysis for more than one diverging terms. Let us take near the 
horizon 

(M 2 = (21) 

where C 2 (x) is a non-diverging function independent of (3 in the leading order over or close to T-L, 
and is independent of R over Ti. 

The divergence of the second term within the square bracket in Eq. (fT8|) comes from {dej) 
which, by Eq. (f2Tj) is C(/3 _1 ). So this term can be neglected with respect to the quadratic term 

(dgJ) 2 . Hence comparing Eq.s (flU)) . (121[) we have D(x) = '"Mf . 

Using Eq. (|15|) we obtain from Eq. ([21]) the following divergence over T-L 

*1 =± 2^ Ci ( X ) (22) 



dRdO, 
On the other hand we can write Eq. (|18[) near % now as 



2 A 2 r 



( X a 9al - eff - D(x) . (23) 



We shall take the Lie derivative of Eq. (|23| with respect to mf over By our choice R a R a = 
A 2 = V a /3 2 V a /3 2 . Also, the function « in Eq. (|14|) is a constant over T-L. This means that c^k = 
over H. Since by our definition the vector field Xh is Killing over TL, the term (xjjda-f — e/) is a 
conserved quantity, i.e., a constant [21]. We shall regard this term to be the conserved effective 
energy (E) of the particle. So, using Eq.s (|T3|) . (f2"2"j) the Lie derivative of Eq. (|23|) with respect to 
m" gives the following 0(/3 _1 ) divergence over H 

d ei D(x) = ±^Ci(x) [E 2 - D(x)] i . (24) 

Eq. (|24|) contradicts the fact that D(x) is independent of (3, A or R in the leading order over H. 
So, Eq. (fT9|) cannot be true. Similarly we can show that the term r ) a bd a Id b I cannot be divergent 
as 0(/r n ) for any n > 2. Thus p 2 j ab d a Id b I = over the horizon. 

With all these, we now integrate Eq. (TTg|) across the horizon along a complex path 

I± = ± / ^(xW-ef)dR, (25) 

where complex integration is understood. The +(— ) sign stands for outgoing (incoming) solution. 
Eq. (|25|) gives the emission (absorption) probability for a stationary black hole satisfying the 
assumptions we have made. 

In order to verify the validity of Eq. (I25|) . at this point we need some particular metrics. We 
shall find out the vector fields Xh an d R a , an d then compute I± from Eq. ([25]) . 
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Let us start with four dimensions by considering the charged Kerr black hole 

, a A-a 2 sin 2 0, ? 2a sin 2 # (r 2 + a 2 - A) (r 2 + a 2 ) 2 - Aa 2 sin 2 . 
ds 2 = dt 2 i LdtdS + ± '- sin 2 

e s e 

+ ^dr 2 + Ed6» 2 , (26) 

where E = r 2 + a 2 cos 2 6*, A = r 2 + a 2 + Q 2 — 2Mr ; a and Q are the parameters specifying 
rotation and charge respectively. A = defines the horizon (rn). The gauge field of this solution 
is A a = [{dt) a - a sin 2 B{dcp) a ] . 

We first define x° — (9t) a — ^-(9^) a , such that Xa(d$) a — everywhere. Near the horizon 

we have XaX a = — 1 ~ ~ ( r 2 +a 2)2 A ^ a2 sin 2 e — ®" ^ ~ ® over ^ aB horizon which implies x a 
becomes null over the horizon and timelike outside it. 

Over the horizon x a becomes, Xh = (dt) a - ^-{rn){d<h) a = (d t ) a + 3 ° 2 (^) a , which is Killing 
and null. Thus we have specified the required vector field x a which becomes null and Killing over 
the horizon. 

Next we need to find out R a and the parameter R along it. Using the expression x a — {dt) a — 
— (dd,) a , we have x a V /3 2 = everywhere. So we can let R a = V a /3 2 . Then using the expressions 
for ft 2 and the metric functions (Eq. ([25])) we have near the horizon 



dR aP P (r 2 +a 2 ) 4 



W = -nr = V a ^V a /3 2 = . „ ■ 2 4 A + 0(A 2 ), (27) 



where the prime denotes derivative with respect to r. Thus we have found out the norm A 2 (= 
V a (3 2 V a f3 2 ) of the vector field R a which becomes null over the horizon. Also. Eq. (|2~T1) gives near 
the horizon 



Thus we have specified the coordinate or the parameter R along R a . Note that Eq. (1281) implies 
that near the horizon, choosing the vector field R a — V a /3 2 means a coordinate transformation 
r — >• R in the metric (|26p. 

The components of the gauge field ^4 a on the horizon are given by A a Xu = — SPT H 9 , and 

A a (d c h) a — t — = — l rH °, s '. n , ; — rr-. The near horizon contribution comes only from the first one. 

( r H + Q COS "J ( r H+ a j 

Substituting the near horizon norms Y a X a = — 1 ~ -tittti, P a P a = A 2 = , 2 ^ S 2 , 4 A' 2 , and 

\t -\-a ) [r -f-ct ) 

dR = {r2+ ^,2 dA into Eq. (J2SJI we have 

/± = ± / (aW - e/) ^-dr, (29) 

where / = -A a x^ = — r iPl" 3 ■ Eq. (J29J) was first obtained in [23j [24] by explicitly solving the 
semiclassical Dirac equation by method of separation of variables. 

iI(r- H ) ± 2 

The emission (absorption) probabilities are given by ~ e k [5]. We shall not go into 

the details of the complexification of the 'path', the choice of contours and explicit evaluation of 
Eq. ([2"9"|. We refer the reader to [23 HI] for this. Explicit evaluation of Eq. (J30J) and the 
emission (Pe) or absorption (Pa) probabilities give the desired temperature of the emission from 
the exponential behaviour of -p^-. The Hawking temperature is found to be Th = ^pp, where kh is 
the surface gravity of the event horizon. 
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After this, we shall consider some examples from higher dimensions. First, we consider non- 
extremal rotating charged black hole solution of five dimensional minimal supergravity with two 
different rotation parameters in the Boyer-Lindquist coordinates [25) . 



ds 2 



Ag (1 + g 2 r 2 ) A 2 , (2mp 2 - q 2 + 2abqg 2 p 2 ^ 



P 4 ^ 



dt 2 



P 



dr 2 



£de 2 

A 9 



(r 2 + a 2 ) sin 2 9 a 2 (2mp 2 - q 2 ) sin 4 9 + 2abqp 2 sin 4 9 



n 4 T 2 
P 



(r 2 + b 2 ) cos 2 9 b 2 (2mp 2 - q 2 ) cos 4 9 + 2abqp 2 cos 4 9 



P 4 ^l 

2A e sin 2 9 [a (2mp 2 - q 2 ) + bqp 2 (l + a 2 g 2 )] 

„4V2 



dip 2 



P 4 S 2 E b 

2A e cos 2 9 [b (2mp 2 - q 2 ) + aqp 2 (f + b 2 g 2 )] 

2 sin 2 6> cos 2 6> [ab (2mp 2 - q 2 ) + qp 2 (a 2 + b 2 )] 
p 4 S a S b 



dtd(j) 



dtdip 
d4>dip, 



(30) 



where p 2 



a 2 cos 2 + b 2 sin 2 0) , A e = (l - a 2 g 2 cos 2 9 - b 2 g 2 sin 2 9) , S a = (1 - a 2 5 2 ), E 6 



(f - &V) and A r = 



(r 2 +a 2 )(r 2 +b 2 )(l+g 2 r 2 )+q 2 +2abq 



2M 



The black hole event horizon (ru) is 



given by A r (rH) = 0. The parameters (M, a, 6, q) specify respectively the mass, angular momenta 
and the charge of the black hole, g is a real positive constant. The gauge field corresponding to 

the charge q is given by A a = ^ (^(dt)„ - ^f^(#)a - . 

The angular velocities of the comoving observers on the horizon are given by [26] 



f2„ 



{g<f><f>9iW> - (ffV</>) 2 } 

{gt-ipg<t><t> - gt<t>94>-4>} 



a{r 2 l + b 2 ){l+g 2 r 2 l ) + bq 
(r| + a 2 )(r 2 I + b 2 ) + abq ' 

6(r 2 +a 2 )(l+ff 2 r 2 ) + ag 
(r| + a 2 )(r 2 I + 6 2 ) + abq ' 



We note that the vector field 

{gt<t>gipijj - gti>g<j>ip} 



x 



(dtT 



{g<t>4>g^ ~ (g^<t>) 2 } 



<q Y _ {gtipg^ - 9t^>g^} , g , a 
21 " {9W9W - (g^<p) 2 } % 



(31) 



(32) 



is orthogonal to (d^Y and (d t f,) a everywhere. Also, the near horizon norm of x a is X a Xa = —ft 2 = 
- i(r 2 +a 2 f(r 2 ?b 2 )+ab q ]' 2 + °( A r)- Thus X a becomes null over the horizon. 

Also, Eq. (f3Tj) shows that x" becomes a Killing field Xh over the horizon, where 



(33) 



So, we have specified the required vector field x a which becomes null and Killing over the horizon. 

Also, exactly through the same manner as in the Kerr-Newman metric, we can specify the other 
null vector field R a , its norm A 2 , and the coordinate R for the metric ([30]). Choosing R a — V a /3 2 
we have near the horizon 



R a R a = A = V a P 2 V a p 2 = R a V a P = 



A r p 2 r 



2 ,A 



\(r 2 + a 2 )(r 2 + b 2 ) + abq]' 



(r 2 A r )' 2 , 



(34) 
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which becomes null over the horizon. Also, near the horizon the coordinate R along R a is given by 



f [(r 2 + a 2 )(r 2 + b 2 ) + ab q ] 2 d(r 2 A r ) 
J (r 2 A r )(r 2 A r )' 2 ' 

The gauge field A a has three components : (A a x^ A a (d c f ) ) a , A a (d^) a ), of which the near horizon 
contribution comes only from A a x%^- 

Substituting the near horizon norms B 2 = „ , Jk \ — nil A 2 = ,, „ , „^. T S , r ,„, , , l4 (r 2 A r )' 2 

[{r 2 +a z ){r 2 +b 2 )+abq\ 1 {{r 2 +a 2 ){r 2 +b 2 )+abq\ V / 

A J T> [(r 2 +a 2 Kr 2 +b 2 )+ab q } 2 d(r 2 A r ) 

and dR = 1 p A r )(r 2 A r )>' 2 mto ^Q- <125J) we have 

/" \(r 2 + a 2 )(r 2 + b 2 ) + abq] 

where / = — A a x^ = — ( r 2 +a 2)(>2^!{,2) +abg ■ Eq. ([36]) was first obtained in [26] by explicit solution 
of the semiclassical Dirac equation by method of separation of variables. Complex integration of 
Eq. ([36]) across the horizon and computation of the emission (absorption) probabilities give the 
expected Hawking temperature in terms of the Killing horizon's surface gravity [25[ I26j . 

It can be easily verified using the same methods as above that Eq. ([25]) also applies well and 
recovers the desired results for the 5 dimensional stationary solutions with Killing horizons like 
Kerr-Godel black hole [27], squashed Kaluza-Klein black hole [28] [29], a black string [28] [30], black 
hole solutions of z = 4 Horava-Lifshitz gravity [31] [32] and the toroidal black hole solutions like in 
[33]. 

Our scheme also applies very easily to an n dimensional generalization of the Kerr black hole 
with a single rotation parameter ^34} 

u 2 r"~ 5 £ 
ds 2 = -dt 2 + (r 2 + a 2 ) sin 2 6d<j) 2 H Hr— (dt - a sin 2 6d(£) H ^—r, ^ dr 2 

+ me 2 +r 2 cos 2 <MT- 4 , (37) 

where the parameters (/x, a) represents the mass and angular momentum of the black hole. S = 
r 2 + a 2 cos 2 9 and dVl n ~ A represents the metric over an (n — 4) sphere. 

Eq. (|25[) applies to a de Sitter horizon also, provided the assumptions stated at the beginning of 
this section are true for that case. Such an example is the Kerr-de Sitter spacetime. The de Sitter 
horizon for this spacetime is a Killing horizon |35j . One can show, following exactly the similar 
way as before that all the other assumptions are valid for this case. Explicit evaluation of Eq. ([25]) 
gives the expected thermal character of the incoming radiation. 

4 Vector, spin-2 and spin-| fields 



Now we shall show that all the approaches and conclusions made in the preceding sections also 

3 
" 2 



hold for the Proca, massive spin-2 and spin-| fields. Let us first consider the equation of motion 



for a Proca field A 



b 



h 

where F a b = V a Ab — VbA a . Eq. ([38]) can be written as 



V a F ab = -—A b , (38) 



V a V<M h - R b a A a - V b {V a A a ) = --oA b . (39) 
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But Eq. (|38j) implies that V a A a — identically. Now let us choose a set of orthonormal basis 
|ei M ' l |. We expand the vector field A a in this basis, Ab = e^A^y With this expansion and the 
fact that V a A a = 0, Eq. j38j becomes 



e^V VM (M) + A (M) V a V Q e ^ + 2V a A (M) V a e^ - R^A W = -^A ( ^\ (40) 
which, after contracting both sides by e^y. , reduces to 

2 

V.VM (l/) + A (fl) e b {v) V a V a e^ + 2e\ v) V a A M V a e^ l) - fl M WA M = ~^\,y (41) 

We choose the usual WKB ansatz for each At v \ : At v ^ = f v {x)e~ ~ , substitute into Eq. (|4"Tj) . and 
take the semiclassical limit ?i — > 0. 

Then it immediately turns out that in the semiclassical limit Eq. (|4ip can be effectively repre- 
sented by n Klein-Gordon equations for the scalars Ai v ^ 

2 

m 

V a VM M + -^A {v) = 0, (42) 

with v — 0, 1, 2, . . . , (n — 1). When each of the Eq.s (|42]) is explicitly expanded and the near 
horizon limit is taken, we get back Eq. ([25]) with e = 0. 

Next, we turn our attention to the massive spin-2 field ir a b satisfying Pauli-Fierz equation |36j 

o 

Tfl 

V c V c 7r a6 + -2-7r a6 = 0, (43) 
n 

where n a b are symmetric tensor fields. As before we expand ir a b in orthonormal basis, ir a b = 
e ^ ^ {p.){v) ■ I n the semiclassical limit and for the WKB ansatz Eq. (|4"3")l can effectively be 
represented by "(" 2 +1 ) Klein-Gordon equations for the scalars TT(p)( v ) 

2 

in 

V c V c tt (m)M + -^n w{v) = 0, (44) 

and thus similar conclusions hold for this case also. 

Finally, we wish to briefly address the spin-| fields satisfying the Rarita-Schwinger equation 
[37]. The tunneling phenomenon for this field was addressed in [35] for the Kerr black hole by 
explicitly solving the equations of motion in the near horizon limit. 

The Rarita-Schwinger equation in a curved spacetime reads 

i7 a V a * 6 = — * b , (45) 
n 

where "ft, = is a spinor with s being the spin index. The 7's are matrices (with matrix 
indices suppressed) satisfying the anti-commutation relation similar to the Dirac 7's: {7°, 7 b } = 
2g ab l. The spin-covariant derivative V is defined as V a ^b := (d a + T a ) ^b, where T a are the spin 
connection matrices (with suppressed matrix indices). Also, ^b satisfies an additional constraint 
7 a *a = 0. 

Due to the similarity of the spin-| fields with the Dirac spinors discussed in Sect. 2, we shall 
apply the same method here to show that ^b satisfies the Klein-Gordon equation in the semiclassical 
WKB framework. So, we square Eq. (|4"5)) by applying i"f c V c from left. A little computation, using 
the definition of the spin-covariant derivative V a , the anti-commutation relation satisfied by the 
7's, and also the commutativity of the partial derivatives yields 

1 m 2 

v„v a * 6 + - [ 7 a , 7 C ] (d [a r c] + r [a r c] ) * fc + ( 7 c v c7 a ) v a * 6 = — -r*&- (46) 



10 



So, as in the previous cases, it immediately follows then for the usual ansatz 



*a = 



A a (x)e— s— 
B a (x)e n 
C a (a;)e-^ 



(47) 



Eq. (|46p reduce to the Klein-Gordon equations in the semiclassical limit. We can easily generalize 
this result for a charged spin-| particle coupled to a gauge field by replacing the spin covariant 
derivative by the gauge spin covariant derivative. This gives charged Klein-Gordon equations. 



5 Discussions 

We now summarize our results. The objective of this work was to put the complex path approach 
for stationary black holes in a general framework. To do this, we have dealt with some well known 
physical matter equations and shown for any arbitrary spacetime in a coordinate independent way 
that in the semiclassical limit the WKB ansatz implies that all those equations of motion are 
equivalent to the Klein-Gordon equation. We have done this without choosing any particular basis 
of the vector fields or the 7 matrices. We needed to assume only that a metric g a b can be defined 
on the spacetime which guarantees the existence of the orthonormal basis | e o^} [21] . So it is 
clear that as far as the semiclassical level is concerned it is sufficient to work only with scalars 
for any arbitrary black hole. It also becomes clear from that the Hawking temperature is indeed 
independent of the particle species we are concerned with. 

We further presented a general coordinate independent expression for the emission probabil- 
ity from an arbitrary stationary black hole with some assumed geometrical properties (Eq. (1250 . 
We showed that finding the emission probability or the Hawking temperature for such black holes 
reduces to merely finding a null coordinate or a null vector field (which is spacelike outside the hori- 
zon), and the norm of the timelike vector field which is orthogonal to the horizon and becomes null 
and Killing over the horizon. At this point we can use any specific metric for explicit computation 
and we illustrated the validity of Eq. (|2~51) by taking several examples. 

The principle message of this work is the following. The semiclassical method provides us a 
way through which we can treat the equations of motions of different spin fields and compute the 
single particle emission probability or the Hawking temperature for a stationary black hole in an 
identical footing or manner. 
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